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THICK RECTANGULAR PLATES ON 
AN ELASTIC FOUNDATION 


Daniel Frederick,* J.M. ASCE 


ABSTRACT 


The governing equations and solutions for the bending of thick rectangular 
plates resting on an elastic foundation are presented using the theory of plates 
developed by E. Reissner. A discussion of various types of fixity of edges in- 
volving three boundary conditions is presented and it is shown that Levy and 

Navier type solutions can be used. Also, two simple examples are included. 


INTRODUCTION 


The thick plate theory developed by E. Reissner’, *, *» in 1945 is significant 
because it includes the effects of shear deformation and normal pressure. In 
addition, it permits the specification of three boundary conditions on an edge 
and requires the evaluation of just two functions. A theory analogous to that 
of Reissner’s was developed by Mindlin‘* in 1951 for the motions of plates. 
Since its development, Reissner’s theory has been applied to the solution 
of several problems in the bending of piates. In 1952, Schafer® presented an 
example of the bending of a rectangular plate under a sinusoidal loading. 
Later, Naghdi and Rowley’ extended the theory to include an elastic foundation 
which behaves according to the classical Winkler-Zimmerman assumption. 
A detailed history of the applications of the Reissner theory is contained in 
reference 7. 
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This paper discusses the thick plate equations and their solutions for 
problems concerning the bending of rectangular plates resting on an elastic 
foundation. The basic equations are those obtained from the Reissner theory 
modified to include an elastic foundation in the manner of Naghdi and Rowley. 
A few simple numerical examples are also presented. 


The Fundamental Equations 


In order to conserve space, Reissner’s notation will be used throughout 
this paper. For a detailed discussion of the notations and assumptions, the 
reader should refer to reference 3. 

Fig. 1 shows the corrdinate system and the stress resultants M,, M_, 
Mxy, Vx, Vy in their positive directions. The displacements u, v and Ware 
positive when they act toward the positive ends of the x, y and z axes respect- 
ively; and they are assumed to be of the form u = z a(x, y), v = zB (x,y) and 
w = wix,y). Here aand fare the rotations in the x and y directions respect- 
ively, of a line element originally perpendicular to the neutral surface and W 
is the vertical displacement of all points on this line element. They are de- 
fined as 


h $ h 
Zz 2 
12 12 4z 
a= uzdz, vzdz, ws w(l - )dz. 
_h 
Zz 2 
Furthermore, it will be assumed that the stresses T xz, Tyz 2nd 07 satis- 


fy the equilibrium equations expressed in terms of the partial derivatives of 
the stress components and with no body forces and are therefore given by 


3V 42° 


In these equations q, is the pressure on the upper face and q, = - kw (k is 
the foundation modulus) is the pressure on the lower face of the plate. The 
stress resultants are defined as in the classical theory and are quantities per 


unit of length. 


Basic Equations Governing the Bending of a Rectangular Plate on an Elastic 
Foundation. 

These equations arise from the application of Reissner’s variational 
principle* taking into account the assumptions mentioned above. They also 
include shear stresses T,,-7 and T,77 on the top face of the plate and the 
shear stresses Tyyp and Ty,, on y - bottom face of the plate. After some 
rearranging, these equations may be written as 


a 2 

(2) 
q q 
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as 


The first six equations are the governing equations which contain six un- 
knowns wW, Vy Vy M,, M, and Myxy- 
Boundary Conditions. 

In addition, the boundary conditions also arise from the variational prin- 
ciple and are three in number for each edge of the plate instead of two as in 
the classical theory. These were given by Reissner', *» * and applications 
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+ 
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S =f S =f - 
x xzT xzB y yxT yo. 


involving a particular edge have been published by several authors. Reissner’ 
has discussed a free edge, Schafer® describes one type of simply supported 
edge and Hirsch® gives boundary conditions at the edge of a rigid circular in- 
clusion. Reference 7 summarizes the applications in which different boundary 
conditions are used. 

Since three boundary conditions are to be satisfied, the usual definitions 
for the various types of fixity of an edge (simply supported, clamped, etc.) 
must be extended. These are discussed below. For the purpose of discussion, 
a straight edge perpendicular to the x axis will be used although the same 
statements apply to an arbitrary curvilinear boundary. The Reissner bound- 
ary conditions are (1) M, = Mx or a@ prescribed, (2) Mxy = = My xy OF B pre- 
scribed and (3) V, = V, or W prescribed. (Primed quantities are those 
prescribed on the boundary. ) 


(a) Simply Supported Edge. 

In addition to the two usual conditions W = 0 and M, = 0, a third condition 
on either B or M,, must be specified. This may be either B = 0 or Myy = 0. 
Since M,_ is not Zero on a simply supported edge in the classical theory, the 
classical’ definition is perhaps more closely approximated by choosing B = 0. 
Physically, this condition may be achieved by using a knife edge on the top 
and bottom of the plate if there is sufficient friction between the knife edges 
and the plate to prevent rotation of the line elements in the y direction. This 
scheme will prevent uplift at the corners which is annulled in the classical 
theory by using concentrated force on the corners. The other possibility, 

= 0, would require that there be no friction between the knife edges and 

the plate. This last set of conditions was used by Schafer® in his example of 
a rectangular plate under a sinusoidal loading. 


(b) Free Edge. 

The free edge was used by Reissner* in an application concerning an in- 
finite plate with a circular hole free of stress. Here, the two conditions of 
the classical theory 


are replaced by M, = 0, Mxy = 0 and Vx = 0. It is no longer necessary to 
transform the twisting moment into an equivalent vertical shear as in the 
classical theory. The vertical shear and the twisting moment each must 
vanish separately. 


(c) Clamped Edge. 


The usual conditions W = 0 and = 0 will be replaced by (1) w = 0, 


(2) a =0, or & a. 0 and (3) B=0 or Mxy = 0. Although the clamped condi- 


tions may sia never realized in an experiment, it appears that the condition 


a@ = 0 is closer to fulfillment rather than the condition ~ =0. Also, since 


the statements made in the paragraph on the simply supported edge regarding 


the third condition, 8 = 0 or Myy = 0 still apply, a clamped edge is perhaps 
best defined by W = 0, @ = 0 and B= 0. 


8. Hirsch. R. A. “The Effect of a Rigid Circular Inclusion on the Bending of 
a Thick Elastic Plate.” Journal of Applied Mechanics, Vol. 19, 1952. 
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(d) Elastically Restrained Edge. 

If a thick elastic plate frames into a beam, for which the governing equa- 
tions are the usual classical ones which do not include the effect of shear 
deformation, then the boundary conditions for the plate can be written as: 


* 


(M, PLATE * 


(V = 


In these equations, B is the flexural and C is the torsional rigidity of the 
beam. These express the conditions that (a) the moment M, in the plate must 
be associated with the twisting moment per unit length in the beam, (b) the 
twisting moment M,, in the plate is equal to the bending moment per unit of 
length in the beam and (c) the vertical shear in the plate must be the load per 
unit of length on the beam. If the theory governing the bending of the beam 
includes shear deformation, then the beam terms would be modified 
accordingly. 

Other types of edges may be reconsidered in the same manner. 


A Method of Solving the Basic Equations.® 
The basic equations (3) will now be considered for the case where 


sy =0. Solutions for this system of equations can be found from 


the equations on W, V,, and Vy. The general solution for W will be written as 
the sum of a particular integral and a complementary function as 
W = Wy + Wp, where 


10 

Similarly, V, and V. can be written as the sum of a particular and com- 
plementary function as Vy = Vxq + Vxp and Vy = VyH + Vyp- If the function 
W (x, y) is introduced, so that 


with the constant c? = 


ay 
Vis By (6) 


then the governing equation for is 


9. Frederick, D.: “On Some Problems in the Bending of Thick Plates.” Ph.D. 
Dissertation, University of Michigan, 1955. 
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The particular solutions for Vx and Vy can be written as 


* a 2— 


V+ 
y By | D +c kw,,) 


where V x and Vy are functions which must satisfy the following equations by 
virtue of the governing equations on V, and Vy: 


22 2 


h 


2 2 2 


To vy =D ty ~ -») p +4). 


Hence, for the rectangular plate, the complementary functions can be 
found from equations (4) and (7) and the constants of integration associated 
with these can be evaluated from the 12 boundary conditions which arise from 
the four edges of the plate. 


Applications 


In this section, two simple examples are presented and the results are 
compared numerically with the predictions of the classical theory. In addi- 
tion, some general solutions are discussed. 


Navier Type Solution.” 

A rectangular plate of dimensions a by b has coordinate axes along two 
edges as shown in Fig. 2. It is simply supported on each of the four edges in 
the following sense: 


M,(0,¥) = M,(a,y) = M,(x,0) = M,(x,b) = 0 
w(0,y) wia,y) W (x, o) = w (x, b) =0 


Blo,y) = 6B (a,y) = a (x,0) = a (x,b) = 0, 


Furthermore, it is assumed that the top and bottom faces are free of 
shearing stress and that the loading function q, (x,y) can be represented by a 
double Fourier Series. 


= sin sin . (11) 


The Fourier coefficients ann @re determined in the usual way 


10. Timoshenko, S.: “Theory of Plates and Shells,” First Edition, McGraw- 


Hill Book Company, 1940. 
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a b 
q, (x, y) sin — sin dxdy. (12) 


With the assumptions above w, V, and V, can be expressed in the Navier 
form and the double Fourier coefficients can be found in the usual way. The 
final expressions for w, V, and Vy are presented below. 


k {,mt 


* 


m=l 


we 


c sin 
y mn 
m=l n= 


The solutions given by equations (13) and (14) satisfy the boundary conditions 
(10) and hence the problem is solved. The classical Navier solution satisfies 
the classical plate equations which are obtained by dropping the terms con- 
taining h? in the governing equations (3). It also satisfies the classical simple 
support boundary conditions which are the conditions (10) with the conditions 
on and omitted. 


Levy Type Solution.*° 

The plate of the previous paragraph will be considered again with a differ- 
ent set of boundary conditions. It will be simply supported on two opposite 
edges in the same manner as above. 


10. See p. 6 ftn, 


© 
w > sin sin (13) | 
m=l n=l 
mn 
a=l (14) 
with 2 2 2 
mt mt 
mn 
+ 
nt mm,” as.” + k + | 
h mt nt 
To } + (>) +1 
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M,(o,y) M, (a, y) =0 
Wlo,y) = Wla,y) =0 (15) 
B(o,y) = Bla,y) = 0. 


Furthermore, it may have any combination of boundary conditions on the other 
two edges consistent with the Reissner theory. In addition, it will be assumed 
that the function q, is independent of y and can be represented by the Fourier 

series 


(x) a, sin, (16) 
n= 


Again, the Fourier coefficients can be evaluated in the usual way and are 


_h nx 
q(x) sin —— dx. (17) 


Due to the even derivatives occuring in equation (5), it is obvious that the 
particular solution Wp can be written in the following form: 


w,= sin (18) 
n * 


Upon substituting equation (18) into equation (5) the coefficients w, are found 
to be 


2 
1 a 
w =- 
n 


DIA) + k 


With Wp, Vx can be determined from the first of equations (9). This can be 
expressed in the same form since the left side of the governing equations con- 
tains the function and even derivatives of the function. 


From the second of equations (9), it can be seen that V* = 0, since the right 
side contains derivatives with respect to y while W Wp arid q, are independent 
of y. 

With the boundary conditions (15) in niind it is evident that for this problem 
Wy can be taken in a Levy form. 
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n= 
where 
1+ 


Upon substituting this into equation (4) and solving the resulting differential 
equation on Y,(y), Wy becomes 


Wy = ye, cosh cos fy + cosh sin fy + Cc. sinh cos 
n=l 


‘ nwx 
+ sinh sin fy) in (21) 


where A,, B,, Cy, and Dy, are independent constants of integration and where 


l l 


nv 
+ 
Again, in view of the boundary conditions (15), W will be taken as 


Y = sin (22) 


n=l 


After substituting this into equation (7) and solving the resulting differen- 
tial equation on Y,), we find that the final form for Wcan be written as 


p (x,y) = 2 | G, cosh B.y +H, sinh By sin“ (23) 


where BY? = a, + eg and G,, and Hy, are independent constants of integration 


With W,, and W given by equations (21) and (23) respectively, the boundary 
conditions (15) are automatically satisfied. There still remains six independ- 
ent constants which will be determined by the six boundary conditions on the 
two remaining edges of the plate. In the classical theory at the corresponding 
stage, the Levy solution has four constants which are determined by four 
boundary conditions on the two remaining edges of the plate. 

The method described above also works for the case where two opposite 
edges of the plate have the boundary conditions 


4 
| 

| 

with | 
k 

| 
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a(o,y) = a(a,y) = 0 
Vfo,y) = Vi =0 (24) 
M fo, y) = M (a, y)=0, 


while the remaining edges have any combination of conditions. In this case, 
the solutions for Wu and Ware taken as products of a function of y and 


cos ae . The boundary conditions (24) apply to an infinite plate supported 


by equidistant columns and this problem also has a counterpart in theclassical 
theory. 


A Line on an Infinite Plate Resting on an Elastic Foundation.*® (See p. 5 ftn.) 

A line load of intensity P (force per unit length) acts along the y axis as 
shown in Fig. 3. Since the variable y does not enter into the solution, the 
governing equations (3) become 


Due to the load, there will be a uniform deflection along the y axis. This 
will be denoted by w,. Hence, the boundary conditions at the origin will be 
written as 


wo) = 
a(o) = 0 (26) 
B(o) = 0 

The general method of solution described above in which a function W was 


introduced will be used. For this particular example since W is independent 
of y, the equation on W becomes 
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2 
dw k= % 4& 
x dx dx 
dV, 
=0 
10 “ax y 
2 dV 
M = pow, h“v kw + 
(25) 
ox 
B= Vy. 


a“ 

- = 0, 27 
y (27) 


Using the condition that w, M, and Vx vanish as x approaches infinity, the 
solutions for YW and w which are to be used take the following form 


x 
cos f x + Be sinf x 


Upon application of the boundary conditions (26), the constants A, B and G 
become 


G =0 


+ {De - 3De - | 


3 
Dit -3Def f 


With these constants, the expressions for M,, V, and can be found from 
equations (25), It can also be seen that M,, and Vy are zero. The above 
solution will now be compared with the classical solution where the boundary 
conditions are 


B=-w 


= dw 
and the governing equations are those obtained from equations (25) by 
dropping the terms containing h*. For convenience in the calculations, it has 


4 
been assumed that = = 1 for both theories. The deflections, moments and 


shears are shown in Fig. 4. In these and the figures which follow C denotes 
the results of the classical theory and R denotes the results of the Reissner 
theory. 

As an alternative, the results may be expressed in terms of P instead of 


Wo through the relation Vg&(o) = .. Fig. 5 shows the deflection, moment, 


shear and the stress g, in terms of the load P. 
Since the bending stresses and shearing stresses in the plate are given by 
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y = Ge 
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e, =F coss, o sing 
1/2 

h 
r = (.-), 6 = tan 
D k 
c 
TD 
A=w 


the usual elementary strength of materials formulas, it can be seen from 
Figs. 4b and 4c that the maximum bending stress is reduced about 55% and 
the maximum shearing stress is lowered by about 15% over the corre- 
sponding values given by the classical theory. Thus, a design of the slab 
based upon the classical theory is “safe” but overdesigned. 

From Figs. 5b and 5d it is seen that the same conclusions are valid for the 
case of a given line load pressing on the Reissner plate and a classical plate. 
The moments and shears are reduced by using a design based upon the 
Reissner theory. On the other hand, a glance at Fig. 5a shows that the maxi- 
mum deflection by the Reissner theory is greater than that given by the 
classical theory. These trends appear to be true in general. In the particular 
case of the maximum deflection, it seems intuitively obvious that a theory 
including shear deformation should give a deflection larger than that given by 
a theory which neglects shear deformation. 

Fig. 5d shows that the normal stress 0, decreases as the distance from 
the load increases. 


Bending of a Plate Due to a Uniform Shear Stress on the Top and Bottom 
Faces with Two Opposite Edges Clamped and of Infinite Length Parallel to 
These Edges. 

To illustrate a boundary value problem involving shear stresses on the top 
and bottom faces, the plate shown in Fig. 6 will be considered. In order to 
simplify the problem as much as possible, the plate has been taken of infinite 
length in the y direction and the shear stress on the top and bottom faces has 
been taken as a constant. In addition, the foundation has been omitted since 
it would change only the form of the solution of W. In practice this type of 
bending might occur in a bin full of material in which there is friction at the 
wall of the bin due to packing of the material. This friction would produce 


both bending and extension of the plate; however, this example illustrates a 
method for handling the bending only. 

Since all quantities involved are independent of y, the governing equations 
as obtained from equations (3) are 


s* =S" =0, s* = 29 
x 
(30) 
2 
x dx 
n2 
— 
dv 
1 - v) x 
dw 
M A 
y dx 
av (30) 
M,, dx 
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dw_ 6 
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For the edges which are clamped, the boundary conditions will be taken as 


W(o) = w(L) = 0 
a(o) = a(L) = 0 (31) 
B(o) = B(L) = 0 


Following a procedure similar to that used in the previous example, the solu- 
tions for W and W can be written as 


2 
= 6.2 


2 3 


y =a, sinh ta, cosh x 


There are now six constants a), a,, a), a3, a,, a4, which are to be determined 
by the six boundary conditions (31). Omitting the detailed work, the final 
values of the costants are 


T 


Ga’ 


Hence, the expressions for w, M, and V, become 


2 2 3 
w=3 (FE) - + 20%) 


v= m= 9 
x Ta’ x 


A plot of these equations for V =} is shown in Fig. 7 for different 2 


values. This problem can not be solved by the classical theory since the 

solution W would be the same as above, and the boundary conditions 

on = dw dw 

wo) = = = . 

of integration be zero. However, by letting L > with Lim (7h) remaining 
h>0 

finite, one can get the solution for a plate of zero thickness subjected to a 

moment of magnitude 7th at every point of the plate from the Reissner theory. 

These results are also drawn in Fig. 7. The bending moment M, for h = 0 is 

zero everywhere. 


= 0 would require that all four constants 


| 
sop Vy: | 
a, = 44 = a5 = 0 
(33) 
2 3 Ga’ 3 3 GLa 
where 
L 12h? | 
| 
(34) 
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CONC LUSION 


The thick plate theory which was developed by E. Reissner, which includes 
the effects of shear deformation and normal pressure, can be used to solve 
certain problems in the bending of rectangular plates resting on an elastic 
foundation. In this paper two simple examples are presented and some results 
are plotted graphically. From the graphical plots, the thickness-length para- 
meter at which the effects of shear deformation and normal pressure become 
important can be determined. 
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Figure 1 - An Element of the Plate for Rectangular Coordinates 
Showing the Stress Resultants. 
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Figure 2 - A Rectangular Plate with Coordinate Axes, 
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Figure 3 - An Infinite Plate with a Line Load of Intensity P. 


y 
Z 
x 
818-16 


Figures 4a, 4b, 4c - A Comparison of the Deflections, Moments 
and Shears in terms of Wo for the Line Load Problem, 
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Figures 5a, 5b, 5c, 5d - A Comparison of the Deflections, 


Moments, shears and _# in Terms of P for the Line Load 
Problem, 
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Figure 6 - A Plate with Shear Stresses on the Top and Bottom 
Faces, 
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Figures 7a, 7b, ic A Comparison of the Deflections, Moments 
and Shears for the Plate Bent by Shear Stresses, 
(v= 1/3; the h/L values are given on 
the curves) 
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